Direct/indirect detection signatures of non-thermally produced dark matter 
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We study direct and indirect detection possibilities of neutralino dark matter produced non- 
thermally by e.g. the decay of long-lived particles, as is easily implemented in the case of anomaly 
or mirage mediation models. In this scenario, large self-annihilation cross sections are required 
to account for the present dark matter abundance, and it leads to significant enhancement of the 
gamma-ray signature from the Galactic Center and the positron flux from the dark matter annihi- 
lation. It is found that GLAST and PAMELA will find the signal or give tight constraints on such 
nonthermal production scenarios of neutralino dark matter. 

PACS numbers: 98.80.Cq, 14.80.Ly, 95.35.+d 



I. INTRODUCTION 

While there are lots of cosmological evidence of the 
dark matter in the universe [l|, |2|, its detailed proper- 
ties remain largely undetermined. Requirements for the 
dark matter particle are the following. (1) It reproduces 
the present matter density of the universe. In terms of 
the density parameter, £l m h 2 ~ 0.11 must be satisfied 
where /i(~ 0-70) is the Hubble parameter in units of 
100 km/s/Mpc [|. (2) It is electrically neutral. (3) It is 
cold, which means that its free-streaming length (Aps) is 
not so long as to seed the structure formation satisfac- 
tory, and this requires Afs < 1 Mpc. 

In fact many candidates of dark matter are proposed in 
the framework of physics beyond the standard model. In 
particular, supersymmetry (SUSY) provides interesting 
candidates. If _R-parity is conserved, the lightest SUSY 
particle (LSP) becomes stable and contributes present 
matter density of the universe. Among SUSY particles, 
the gravitino and (lightest) neutralino are possible can- 
didates of the LSP dark matter. From the viewpoint 
of detection possibility, the gravitino dark matter is un- 
desirable because its interaction strength with ordinary 
matter is Planck-suppressed. 1 In the following our focus 
is the neutralino dark matter, which may have distinct 
signatures of direct and/or indirect detection. 

Usually neutralinos are assumed to be produced ther- 
mally as in the following scenario [J, Q. In the early 
universe with temperature T > 1 TcV, SUSY particles 
including neutralinos are thermalized and their number 
density is given by ~ T 3 . As the temperature decreases, 
their thermal abundance receives Boltzmann suppression 
factor and eventually they decouple from thermal bath 
at the freeze-out temperature Tf ~ m^/20 where tol 
denotes the LSP mass. After that the number density 
of the LSP per comoving volume remains constant until 
now and hence contributes as dark matter of the uni- 



verse. The resultant abundance of the LSP is estimated 
as 
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1 Recently it is pointed out that the detection of inflationary grav- 
itational wave background can help the situation |4(. 



where (cry) denotes the annihilation cross section of the 
LSP and Mp is the reduced Planck scale (= 2.4 x 
10 18 GeV). 

However, such a thermal relic scenario does not always 
hold in realistic SUSY models. For example, there often 
exists Polonyi or moduli field in order to break SUSY and 
gives rise to correct order of gaugino masses. Those sin- 
glet scalar fields generally have long lifetime and decay 
after freeze-out of the LSP, yielding substantial amount 
of LSPs. Actually Polonyi/moduli dominate the uni- 
verse before they decay, and hence reheat the universe 
again with very low reheating temperature of O(l) MeV- 
0(1) GeV depending on their masses [g, 0|- In this case 
a large amount of LSPs are produced non-thermally by 
the Polonyi/modulus decay, and hence large annihilation 
cross section is needed to account for the present dark 
matter abundance. Therefore, taking into account non- 
thermal production mechanism may significantly change 
the properties of the LSP and its direct/indirect detec- 
tion signatures. 

Thus in this paper we study direct/indirect detection 
signatures of non-thermally produced neutralino dark 
matter. As for direct detection, there are some ongo- 
ing and planned projects devoted to detect scattering 
signals of the LSP with nucleons such as CDMS @] and 
XENON 0]. As for indirect detection, many possible 
ways are proposed. First, neutralinos accumulated in the 
Galactic Center annihilate each other and produce line 
and continuum gamma rays. Such gamma-ray signals 
can be searched by satellite experiments (GLAST |10 |) 
or g round-based Cerenkov telescope (HESS [ll[, MAGIC 
[if, CTA [ljj). 

Second, anti-matter such as positrons or anti-protons 
are produced by the annihilation of the neutralinos. Since 
these particles are diffused by galactic magnetic fields 
during their propagation to the Earth, we need to solve 
its propagation in a diffusion model to discuss their flux 



on the Earth [Tj]. Fortunately, the positron flux is less 
sensitive to the precise diffusion model since magnetic 
fields easily dissipates their energy through the propaga- 
tion and positrons come only from near the Earth. These 
anti-matter signals can be detected PAMELA [15| and 
AMS-02 Q. 

Third, neutralinos trapped in the Sun annihilate and 
produce high-energy neutrinos. Super Kamiokande I17l . 
AMANDA [jj], IceCube [H and planned KM3NeT [Si 
experiments search high energy muon signals, which arise 
from high-energy neutrino interaction with Earth matter. 
We investigate characteristic signals of nonthermal neu- 
tralino dark matter on these experiments. 

For the sake of concreteness, we stick to two SUSY 
breaking models : minimal anomaly-mediated SUSY 
breaking model [21| and mirage- mediation model [22j ]. 
The former model predicts the wino-like neutralino LSP 
in broad parameter regions. Since the wino-like neu- 
tralino with a mass of 0(100) GeV has too large annihi- 
lation cross sections, its thermal abundance becomes too 
small to account for the present dark matter abundance. 
Hence we need to consider some non-thermal production 
processes of the neutralino dark matter. The latter model 
contains a heavy modulus field and non-thermal produc- 
tion of the neutralino dark matter is expected naturally. 
As we will see later, the large annihilation cross section 
of the neutralino dark matter is a general feature of non- 
thermal production scenario, and hence our results are 
less sensitive to the model construction. 

A similar subject was studied in Ref. [23||, where it was 
pointed out that the large annihilation cross section of 
the neutralino yields enhancement of the anti-matter sig- 
nals. We emphasize that such an enhancement is rather 
generic feature when considering the nonthermal produc- 
tion scenario of the dark matter, and its detection may be 
directly related to the early Universe cosmology, in par- 
ticular, the existence of late-decaying particles and their 
decay temperature. Also we have performed more de- 
tailed parameter analyses both in the anomaly-mediated 
SUSY breaking and mirage-mediation models, including 
the gamma-ray signature as well as anti-matter searches. 

This paper is organized as follows. In Sec. |TT] we re- 
view some non-thermal production mechanisms of the 
neutralino dark matter. The minimal anomaly-mediated 
SUSY breaking model and the mirage-mediation mode 
are taken as examples. In Sec. Mil detection possibilities 
of nonthermal dark matter are discussed. These include 
direct detection using recoil of nuclei by the neutralino, 
gamma-ray flux from the neutralino annihilation at the 
Galactic Center, positron flux from the annihilation near 
the Earth, and high energy neutrino flux from the anni- 
hilation in the Sun. Sec. HVl is devoted to our conclusions. 

For calculating these direct/indirect detection rates, 
we have utilized DarkSUSY code fUl- 



II. NONTHERMAL PRODUCTION OF 
NEUTRALINO DARK MATTER 

The LSP neutralino is often referred to as a promising 
candidate of the dark matter of the universe. Although 
the standard thermal relic scenario is often assumed, the 
production processes of the neutralino are not limited to 
it in general. In particular, non-thermal production pro- 
cesses may be relevant for estimating the present dark 
matter abundance. Actually we often encounter the cos- 
mological scenarios which include long-lived matter. If 
the long-lived matter has non-negligible fraction of the 
total energy density at the time of its decay, neutralinos 
emitted by its decay processes may amount to significant 
contribution to the dark matter abundance 25]. Here we 
present examples of such late-decaying matter. 

(1) Gravitino : The gravitino is the superpartner of 
the graviton and its interaction strength is suppressed by 
the Planck scale or SUSY breaking scale. Gravitinos are 
efficiently produced in the early universe thermally [26| 
or non-thermally by the NLSP [23, [28| and/or inflaton 
decay 29] . Its thermal abundance is given by 
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for a heavy gravitino where m 3 1 2 is the gravitino mass 
and Tr is the reheating temperature of the universe de- 
fined as T R = (10/7r 2 ^) 1 / 4 Vr i „fA/p with the total de- 
cay rate of the inflaton r; n f . If the gravitino is unstable, 
it decays at later epoch after the freezeout of the LSP. 
In the case of heavy gravitino (m 3 / 2 > 100 TeV) as is 
realized in anomaly-mediated SUSY breaking, it decays 
well before BBN begins without affecting the success of 
BBN, and produces LSP non-thermally. The lifetime of 
the gravitino is estimated as 
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if the gravitino is the heaviest among MSSM particles. 
Depending on the reheating temperature, the abundance 
of non-thermally produced LSPs by the gravitino decay 
can exceed that of thermally produced ones. 

(2) Polonyi field : A Polonyi field is a singlet scalar 
responsible for SUSY breaking. A Polonyi field is always 
required in gravity-mediation models in order to gener- 
ate sizable gaugino masses. However, generically Polonyi 
has the mass (m x ) of order of the gravitino (m 3 / 2 ) and 
its interaction is Planck-suppressed. Thus decay of the 



coherent oscillation of the Polonyi causes cosmological 
disaster, unless Polonyi decays well before BBN or it is 
diluted by some additional entropy production processes. 
The Polonyi field begins to oscillate when the Hubble pa- 
rameter become equal to the Polonyi mass, and its abun- 
dance is estimated as 
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where \o is the initial amplitude of the Polonyi field. 
Here 7 is defined as 
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where T osc = (90/7r 2 g* ) 1 / 4 ^/m x Mp. If the Polonyi mass 
is larger than ~ 100 TeV, it decays before BBN begins. 
The lifetime is estimated as 



t 1 m£ 



4.9 x 10" 2 sec 



100 TeV V 
m x J 



(6) 



In general, Polonyi field also decays into SUSY particles 
and their abundance may be bigger than the thermal relic 
one 0, 0, [3CJ, [31j . Thus we must take into account the 
abundance of the LSP arising from the Polonyi decay. 

(3) Moduli : Modulus field is a scalar field appear- 
ing in the low energy effective theory of string theory 
when the extra dimensions are compactified. Proper- 
ties and their cosmological effects are similar to those 
of the Polonyi 32]. In some particular models the mod- 
ulus mass is related to the gravitino mass. For exam- 
ple, in KKLT setup described below, m x ~ 47r 2 m 3 / 2 is 
obtained. In such a case the decay temperature of the 
modulus can be as large as ~ 0(1) GeV :33], although 
gravitinos produced by the modulus decay may cause an- 
other cosmological difficulty [34,l35j. 2 

(4) Saxion : In SUSY extension of the axion models, 
there exists an additional light scalar degree of freedom, 
which obtains a mass from SUSY breaking effect, called 
saxion (a) [3 71 ] . The saxion mass is naturally expected to 
be ~ m 3 12 and its interaction strength is suppressed by 
the Peccei-Quinn scale /pq(~ 10 10 ~ 12 GeV). Saxions are 
produced in the early universe via the coherent oscillation 
and their decay process affects cosmology 38, 39, 40] . Its 



2 See also [36i | for models of heavy moduli and their cosmological 
effects. 



coherent oscillation contribution is given by 
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where cto is initial amplitude of the saxion and 7 is defined 
analogously to the Polonyi case. If we assume that the 
main decay mode is a — > 2g, where g denotes the gluon, 
the lifetime is estimated as 
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If the saxion is heavier than the LSP, decay modes into 
LSPs or SUSY particles become open and they give sig- 
nificant fraction of the relic LSP density [4l|. Axino, 
which is the fermionic superpartner of the axion, may 
also produce large amount of LSPs by its decay, if unsta- 
ble gim. 

(5) Q-Ball : Q-ball is a non-topological soliton whose 
stability is ensured by a global U(l) symmetry [44j ]. 
In SUSY, flat direction condensates (called Affieck-Dine 
field) can develop to large field value during inflation and 
coherent motion of the Affleck-Dine (AD) field can create 
observed amount of baryon asymmetry 45] . Through the 
dynamics of the AD field, fluctuation of the AD field de- 
velops if its potential is flatter than the quadratic one and 
then it fragments into solitonic objects, Q-balls [46l.l47j]. 
Here the global U(l) charge required to stabilize the Q- 
ball configuration is baryon number. Assuming that the 
AD field begins to oscillate at H ~ m^, with amplitude 
4>i, the Q-ball charge (Q) is estimated as 



= 7^ 
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where 7~6xl0 _3 isa numerical constant and e is called 
ellipticity parameter which is smaller than one. In the 
case of AD baryogenesis using flat directions lifted by 
non-renormalizable superpotential Wnr ~ <f> e /M 3 with 
cutoff scale M, we obtain fy ~ (m^M 3 ) 1 / 4 . Thus the 
charge of Q-ball becomes 
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Here we assume gravity-mediation type Q-ball. 



Once a Q-ball is formed, its lifetime is determined by the 
charge of Q-ball J49| , 
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where m^ denotes the mass of the AD field and k is a 
model dependent constant which can take the value ~ 1- 
10~ 4 . On the other hand, the abundance of Q-ball is 
given by 

PR. = It, f A" 2 
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Thus decay of Q-balls yields significant amount of 
LSPs non-thermally [50, [51J], and they can even domi- 
nate the universe before the decay [52| |. Interestingly, 
this model can explain the present dark matter abun- 
dance and baryon asymmetry simultaneously by the Q- 
ball decay. 

All the above models predict decay temperature of or- 
der 0.1-1 GeV typically, which is smaller than the freeze- 
out temperature of LSP. In this case the non-thermally 
produced LSP abundance may exceed the thermal relic 
one. Thus it is important to reconsider the abundance of 
the LSP produced by the decay of those long-lived mat- 
ters. Hereafter we collectively denote such a long-lived 
matter as x- The following arguments do not depend on 
the detailed properties of x once the decay temperature 
of x is fixed. 

Now let us write down the Boltzmann equations which 
govern the evolution of the number density of the LSP, 
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(14) 

AHp r = (m x — 2Bjjni)T x n x + mj J (crv)nj j , (15) 



where 7Jl and n x denotes the number density of the LSP 
and late-decaying particle; (av) denotes thermally aver- 
aged annihilation cross section of the LSP 4 ; p r denotes 
the radiation energy density; TOl denotes the LSP mass; 
.Bl denotes the branching fraction of the x decay into 
LSP or SUSY particles; and H denotes the Hubble pa- 
rameter. This set of equations is simplified for t > 1/L X 
after x decays, 



hh + 3JJriL 
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This equation can easily be solved and the resulting 
abundance of LSP is simply given by 
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4 Neutralinos are expected to soon reach kinetic equilibrium due 
to interactions with background particles in thermal bath 0, l5.'j| . 



where we have defined number-to-entropy ratio as Yl = 
nh/s with entropy density s. The initial abundance af- 
ter x-decay Yl(T x ) under sudden decay approximation is 
given by 



Y^(T X ) = 
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where B(x — > 2LSP) is the branching ratio of x m to 

(TP) 

LSPs and F L denotes the contribution from thermal 
freezeout, taking into account the dilution from ^-decay. 
If the annihilation cross section is sufficiently small or 
the initial abundance of the LSP is negligible, LSPs can- 
not annihilate each other after \ decays. However, if 
the annihilation cross section is large enough, the LSP 
abundance becomes inversely proportional to the anni- 
hilation cross section similar to the case of thermal relic 
abundance, 
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The crucial difference is that larger annihilation cross 
section is required in order to account for the present 
dark matter density, because the decay temperature of % 
is smaller than the typical freeze-out temperature of the 
LSP, T f - to l /20. 

In Fig. Q] the resulting LSP abundance is shown for 
tol = 300 GeV as a function of annihilation ross sec- 
tion. Dashed (solid) lines correspond to Y^(Tx) = 
10~ 9 (10~ n ) and left (right) ones correspond to Tx = 
1(0.1) GeV. It is seen that for large annihilation cross sec- 
tion the result becomes independent of the initial abun- 
dance, since the LSP abundance is saturated due to the 
annihilation effects. 

Thus we are interested in models which predict large 
annihilation cross section. For concreteness, we re- 
strict ourselves to two SUSY breaking models : mini- 
mal anomaly-mediated SUSY breaking model (MAMSB) 
and mirage- mediation model (mMSB) . The former model 
predicts the wino-like neutralino LSP, which has natu- 
rally large annihilation cross section compared with the 
case of bino-like neutralino. In the latter model, the an- 
nihilation cross section of bino-like neutralino LSP could 
be enhanced through S'-channel resonance, as we will see. 
Now let us briefly introduce these models. 



A. Minimal anomaly-mediation model 

In the anomaly-mediated SUSY breaking model, the 
effect of SUSY breaking in the hidden sector sequestered 
from the observable sector is transmitted to the MSSM 
sector through super- Weyl anomaly effect. Since the 
SUSY breaking effect is suppressed by a loop factor, 




cross section [GeV" ] 



FIG. 1: Nonthermally produced LSP abundance (Qhsph 2 ) 
as a fuction of the annihilation cross section ((er«)). Dashed 
(solid) lines correspond to Y L (T X ) = 10~ 9 (10~ n ). Left 
(right) ones correspond to T\ — 1(0.1) GeV. Here we have 
taken m L = 300 GeV. 



the gravitino becomes a few magnitude heavier than the 
SUSY particles. For example, gaugino masses are given 



by 
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where g^s are gauge coupling constants for i = 1-3 cor- 
responding to U(l)y, SU(2) and SU(3) gauge groups. 
The beta function coefficients 6,'s are calculated as b\ = 
33/5,62 = 1 and 63 = —3. Note that the expression (|20|) 
is valid for all energy scale /i, not for only the input scale. 
Thus it is seen that for weak scale gaugino masses, the 
gravitino mass can be as heavy as 0(100) TeV. One prob- 
lem of such a scenario is that it predicts tachyonic slepton 
masses. Thus in the minimal anomaly-mediation model, 
an additional universal scalar mass term nig is introduced 
to give positive contribution to the scalar masses, as 
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where (3 g and 0f are the beta functions of corresponding 
gauge and Yukawa coupling, and 7 = d In Z/d In /i with 
Z denoting the wave function renormalization. Therefore 
this model is characterized by the following parameter 
set, 



m 3 / 2 ,mo,tan/3,sign^, 



(22) 



where tan/3 = (H u ) / (Hd) represents the ratio of VEVs of 
up- type and down- type Higgses. Here and hereafter we 
assume /1 > because positive /x is favored from muon 
g — 2 experiments. 

In this framework the wino-like neutralino (W°) nat- 
urally becomes the LSP, as is easily seen from Eq. (|2"0"|) . 



The annihilation process through J-F°VF° — > W + W" is 
not helicity-suppressed and hence it has large annihila- 
tion cross section. The annihilation cross section of this 
process is calculated as [3(| 
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where «2 is the SU(2) gauge coupling constant and mw 
denotes the W-boson mass. Actually the wino mass 
should be ~ 3 TeV to account for the dark matter den- 
sity if it is p roduced in the standard thermal freeze-out 
scenario [5J|, and such a heavy LSP mass seems to be 
disfavored from the viewpoint of naturalness. However 
if we extend the production mechanism of the LSP to 
non-thermal origin, light wino mass of 0(100) GeV is 
favored. 

In Fig. [21 we plot the reheating temperature of the 
X-decay, T x , in the (m 3 / 2 ,"io) plane with tan/3 = 10, 
assuming the present dark matter is produced by the de- 
cay of x held. In the most parameter regions, wino-like 
neutralino LSP is realized with the mass of ~ TO3/2 /400. 
For larger mo, the higgsino mass becomes smaller and 
the LSP neutralino contains much Higgsino components. 
Thus, too large values of mo are excluded by the absence 
of electroweak symmetry breaking, and the Higgsino- 
like LSP is realized near the boundary. As can be seen 
from this figure, we expect non-thermal production with 
0.1 GeV < T x < 10 GeV for the neutralino dark matter 
with the mass of O(100)GeV. 



B. Mirage- mediation model 

The mirage-mediation models are based on recent de- 
velopments on moduli-stabilization mechanism in string 
theory, that is, KKLT construction [55j. We denote 
the modulus as T and assume the following type of 
Kahlcr potential, superpotential and gauge kinetic func- 
tion (hereafter we set Mp — 1 unless explicitly written), 
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where $j denotes MSSM superfields, a and k are real 
constant, and Z,{T + T*) = 1/(T + T*) n \ The scalar 
potential for the modulus is given by 



V T 
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where D T W = W T + K T W and K Tf = (K Tf )~ 1 with 
subscript T denoting field derivative with it. An analysis 
shows that this potential has a supersymmctric anti-de 
Sitter (AdS) minimum. Thus in order to obtain de Sitter 
(dS) vacuum consistent with current cosmological obser- 
vations, the additional uplifting potential is needed. This 
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FIG. 2: Contours of T x in order to reproduce correct 
dark matter abundance and asi for the minimal anomaly- 
mediation model. Red thick solid lines represent T x — 
0.1, 1, 10 GeV from left to right. Dot-dashed lines show the 
lightest Higgs mass nth = 115 GeV and 120 GeV and dashed 
lines correspond to fi — 1, 1.5, 2 TeV from left to right re- 
spectively. Upper left region is excluded from electro-weak 
symmetry breaking (EWSB) constraint and lower right region 
predicts stau LSP or tachyonic slepton masses, and hence is 
excluded. 



is provided by, for example, adding an extra brane which 
is sequestered from the observable brane. The additional 
term is 
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where to is 0(1) constant. Thus total scalar potential is 
given by V — Vr + Vim- After fine-tuning the value of 
D, the desired dS minimum is obtained. 

Phcnomcnologically this model provides a character- 
istic pattern of SUSY breaking effect. SUSY is domi- 
nantly broken by the uplifting term introduced to make 
the vacuum energy positive, which becomes a source of 
anonraly- mediation effect. On the other hand the mod- 
ulus T has non- vanishing F-term (F T ) at the resulting 
vacuunr, which becomes a source of modulus-mediation 
effect. Thus a mixture of anomaly- and modulus medi- 
ation is realized, called mixed modulus- anomaly media- 
tion [5^]. Interestingly, these two effects are comparable 
in general, that means F T /(T + T*) ~ m 3/2 /(47r 2 ). For 
example, gaugino masses at GUT scale are given by 
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where M denotes the modulus-mediation contribution 
to the gaugino masses at GUT scale and calculated as 
Mq = F T cVhiRe(/ a ). Taking into account the one-loop 



renormalization group evolution, gaugino masses at the 
scale n is calculated as 
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where the mirage scale M m ; r is defined as 
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Here we have defined a parameter a which characterizes 
the ratio of the anomaly- to modulus-mediation contri- 
bution as 
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Thus all gaugino masses seem to be unified at the mirage 
scale, which leads to the term of 'mirage'-mediation. In 
the original KKLT setup a = 1 is predicted, but a can be 
regarded as a free parameter in nrore general setup. For 
a = 1 , the intermediate scale mirage unification (M m i r = 
3 x 10 9 GeV) is realized. The case of a = 2 is called TeV- 
scale mirage mediation, since M m ; r ~1 TeV is predicted. 
Tachyonic slepton mass problem in the pure anomaly- 
mediation model is naturally solved in this framework, 
because of the modulus mediation contribution. 

Denoting the modulus mediation contribution to the 
sfermion masses and j4-terms at GUT scale as mf and 
Ai, this model is classified by the following parameters, 



Mq, Ci,a>i, tan/3, a, 



(33) 



where we have defined a = mf/M^ and at = Ai/M , 
which are related to m as c, = ai = 1 — rii . 

As is obvious from the construction, this model pre- 
dicts a modulus field T whose mass is estimated as 
tot ~ (87t 2 )to 3 / 2 . The modulus field is likely to dom- 
inate the universe due to its large energy density stored 
in the form of scalar condensates with a large initial am- 
plitude of order the Planck scale. Thus its decay leads to 
the non-thermal production of the neutralino dark mat- 
ter with a reheating temperature of 



T T ~ 170 Me 



VV~c( 



tut 



10 3 TeV/ 



3/2 



(34) 



where we have used the decay rate of the modulus, 
Ft = cto t /47tMJ, with 0(1) constant c. In principle, 
other non-thermal production processes could occur after 
the modulus decay, and also the non-thernral production 
with low reheating temperature might be circunrvented 
by the small initial amplitude or heavier nrodulus mass. 
Thus, we treat the final reheating temperature T x as a 
free parameter, and pay a special attention for the case 
of T x < Tt which should be satisfied once the modulus 
field dominate the universe in the minimal KKLT setup. 
Fig. [3] shows T x in the (M ,tan/3) plane with a = 1, 
cm = o-m = 1 for matter fields and ch = o-h = for 
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FIG. 3: Contours of T x = 10 GeV and 1 GeV (red thick 
solid lines) for the mirage-mediation model. The middle gray 
shaded region corresponds to T x < Tt GeV. The upper left 
blue shaded region is excluded by b — * S7 constraint. Dot- 
ted line shows mj = 115 GeV and dashed lines represent 
(j, = 300, 500, 700 GeV from left to right respectively. 



Higgs fields in the non-thermal production scenario by 
the x-decay. In these parameter set, LSP is bino-likc 
neutralino with a mass of Mi ~ (0.4 + 0.3a)Mo and its 
annihilation cross section is enhanced by the s-channel 
Higgs resonance for 10 < tan/3 < 20, and these parame- 
ter regions are compatible with the non-thermal produc- 
tion scenario with T x < lGeV. 



III. DETECTION OF NONTHERMAL DARK 
MATTER 

A. Direct detection 

Direct detection experiments attempt to observe the 
recoil energy of target nucleus scattered elastically by 
the dark matter. The scattering of nucleus is discussed 
in two classes of interactions: spin-dependent (SD) and 
spin-independent (SI) interactions. For the neutralino 
dark matter, scattering through spin-independent inter- 
actions tends to become important and promising to de- 
tect the signals. This scattering is mainly occured by the 
exchange of higgs fields and squarks, and the cross sec- 
tion is generally irrelevant to the self-annihilation cross 
section. Hence, the non-thermal dark matter, which have 
large self-annihilation cross section, does not necessarily 
lead to the enhancement of direct detection signals. 

Figs. S] and [5] show the elastic LSP-nucleus SI cross 
section in anomaly and mirage mediation models re- 
spectively. Input parameters are choosen randomly 
with 10 TcV < to 3/2 < 300 TeV, 0.5 TeV < 



too < 10 TeV and 3 < tan/? < 50 for MAMSB, 
and 200 GeV < m 1/2 < 600 GeV, 0.5 < a < 
2, 3 < tan/3 < 50 and {cm-,0'MiC-h,0'h) = 
(1, 1, 1, 1), (1, 1, 0, 0), (0.5, 0.5, 0.5, 0.5), (0.5, 0.5, 0, 0) for 
mMSB. In each parameter sets, several phenomenolog- 
ical constraints are imposed, such as b — ► sj constraint, 
Higgs mass bounds, correct electroweak symmetry break- 
ing condition and neutralino LSP condition, and the cor- 
responding cross sections are scatter-plotted as a func- 
tion of the LSP mass only for phenomenologically viable 
parameters. We assume non-thermal dark matter pro- 
duction scenario and calculate the appropriate reheat- 
ing temperature that explain the present dark matter 
abundance for each parameters. In these figures, red 
and green points correspond to lGeV < T x < lOGeV 
and T x < lGeV respectively. For the mMSB case, blue 
points represent the case that the reheating temperature 
is lower than the modulous decay temperature. 
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FIG. 4: Expected SI neutralino-nucleus cross section in 
MAMSB as a function of the LSP mass. Red and green points 
correspond to 1 GeV < T x < 10 GeV and T x < 1 GeV re- 
spectively. The upper dashed line shows the current experi- 
mental bound by CDMS, and the lower line represents future 
expected sensitivity of SuperCDMS (stage C). 



In the MAMSB model, LSP is wino-like neutralino in 
almost all parameter region, and its annihilation cross 
section is directly correlated with the LSP mass. For 
example, we need T x ~ lGeV for tolsp = 300GeV to ex- 
plain the present dark matter abundance by non-thermal 
production. On the other hand, the LSP-nucleus scatter- 
ing amplitude highly depends on the fraction of Higgsino 
components in the LSP neutralino, which is increased 
for larger too. Hence, the low reheating temperature and 
LSP-nucleus SI cross sectionin are completely irrelevant 
as can be seen in Fig. 2] 

The same argument can be applied to the mMSB mod- 
els. In this case, large annihilation cross sections required 
for the non-thermal dark matter scenario can be real- 
ized through s-channel resonance for the bino-like neu- 
tralino or by the large SU(2) gauge interaction for the 
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FIG. 5: Expected SI neutralino-nucleus cross section in 
mMSB as a function of the LSP mass. Red, green and blue 
points correspond to 1 GeV < T x < 10 GeV,T T < T x < 
1 GeV and T x < Tt respectively. Dashed lines are same as 
those in Fig. 21 



Higgsino-like neutralino. However, the annihilation pro- 
cesses are not directly related to the scattering amplitude 
in both cases, and the predicted LSP-nucleus SI cross sec- 
tion spread in a wide range. In the Fig. [51 Higgsino-like 
neutralinos appears with large SI cross sections around 
ctsi — 10~( 42 ~ 43 )cm 2 . This is because mMSB models ex- 
hibit rather compressed mass spectrum and the Higgsino- 
like neutralino has much bino component. 

Present strongest upperbounds are given by the 
CDMS experiment [57| for dark matters with mass of 
O(100)GeV and the LSP-nucleus SI cross section must 
satisfy ctsi ^ 10~ 43 cm 2 . Thus, Higgsino-like neutralino 
produced at T x < lGeV is disfavored by this bound in the 
mMSB model. On the contrary, wide parameter regions 
are allowed both for the bino- like neutralino in the mMSB 
and the wino-like neutralino in the MAMSB. They might 
be explored around ctsi ^ 10~ 46 cm 2 in the projected fu- 
ture sensitivity of SuperCDMS, stage C. Thus, the possi- 
bility to detect the non-thermally produced dark matter 
by direct detection experiments is highly dependent on 
the model parameters irrelevant to its production pro- 



B. Indirect detection 

1. Gamma-ray flux from the Galactic Center 

Observations suggest that the total mass of galaxies 
are dominated by dark matter. Dark matter forms halo 
around the galaxy although its density profile still has 
some uncertainty. We parametrize density profile of the 



Galactic halo as 

p(r) = 



Po 



fe)1 1+ fe)' 



(35) 



with r corresponding to the distance from the Galactic 
Center. Navarro-Frenk- White (NFW) profile [58J corre- 
sponds toa=l,6=l,c = 3 and isothermal profile cor- 
responds to a — 0, b = 2, c = 2. If a > the profile shows 
cuspy structure at the Galactic Center and hence the an- 
nihilation rate is expected to be enhanced compared to 
cored profile. 

Neutralino annihilation processes produce both 
monochromatic and continuum photons, and there 
are many studies related to this issue [59|, |6||. The 
former originates from, e.g., \X ~^ 77 > %1 but we have 
found that branching ratio into these modes are small 
for interesting parameter regions. Thus hereafter we 
concentrate on continuum gamma-ray flux coming from 
cascade decays of annihilation products, mainly from 
pion decays. 

Gamma-ray flux produced by the neutralino annihila- 
tion at the Galactic Center is expressed as [591 ] 



$y(ll>,E) 



dN-y 



87rm 2 dE 



P 2 (l)dl(^), (36) 



l.o.s. 



where the integration is carried out over the line of sight 
and dNj/dE represents the differential number of pho- 
tons with the energy E produced by the neutralino an- 
nihilation. This expression shows that the density pro- 
file dependent part and particle physics model dependent 
part can be separated out. 
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FIG. 6: Expected gamma-ray flux from the Galactic Cen- 
ter above the Energy of 1 GeV with the GLAST satellite as a 
function of the LSP mass in the MAMSB model. The isother- 
mal density profile is assumed. 



It is convenient to define the dimensionless quantity 
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FIG. 7: Expected gamma-ray flux from the Galactic Center 
above the Energy of 1 GeV with the GLAST satellite as a 
function of the LSP mass in the mMSB model. The isother- 
mal density profile is assumed. 



the gamma ray signal and the non-thermal production 
of the dark matter. Thus, large annihilation cross sec- 
tion required for non-thermal dark matter directly leads 
to enhancement of the gamma ray signal, and its detec- 
tion possibility becomes increased compared to the case 
of thermally produced dark matter. 

GLAST can detect dark matter annihilation signals if 
%{E > 1 GeV) > 10~ 10 cm-V 1 . Thus, in both cases, 
it seems difficult to find the gamma-ray signals for the 
cored profile such as the isothermal one. However, the 
flux becomes a few orders of magnitude larger for more 
cuspy profile, such as the NFW profile, and in that case 
GLAST may detect dark matter annihilation signals. To 
summarize, although the definite prediction is impossible 
due to large uncertainity of the density profile, there is 
large possibility to detect the signals from the neutralino 
dark matter produced non-thermally. 



2. Positron flux 



J(tp) as 



1 



1 



.5 kpc V0.3 GeV/cm 3 



p 2 (l)dl(rP), 



and its averaged value over solid angle Aft, 

1 



{J)A " = AnJ An dnJW > 



(37) 
(38) 



where Ail = 2tt(1 — cos(^ max ))- In fact the actual obser- 
vations have some angular resolution AQ and hence the 
gamma-ray flux from the Galactic Center should be inte- 
grated over the solid angle. Performing angular integral, 
Eq. (13"6")) can be rewritten as 



-J J ■J.,-l°"^7 



$ 7 (£) -2.8 x 10~ 12 cm^s 



dE 



I TeVV 



3 x 10- 26 cm 3 /s 



(39) 



(J) A n A^ 



Besides the annihilation cross section and the neutralino 
mass, the gamma-ray flux crucially depends on the di- 
mensionless quantity (J)& a A£l, which is solely deter- 
mined by the density profile of the dark matter halo. 
This is numerically calculated for each profile and solid 
angle. For the GLAST, (J) An AQ = 3 x I0~ 4 and 0.1 for 
the isothermal and NFW profiles respectively. 

In Figs. [5] and [7] expected integrated 7-ray fluxes from 
the Galactic Center $ 7 (-E > I GeV) are shown for 
MAMSB and mMSB respectively. Here the isother- 
mal density profile is applied conservatively. Green and 
red points correspond to lGeV < T x < fOGeV and 
T x < lGeV respectively, and blue points is T x < Tt- 
As is apparent from the figures, the gamma ray flux is 
inversely proportional to the reheating temperature be- 
cause the same annihilation process is relevant both for 



Similar to the gamma-rays described above, positrons 
are also yielded by the neutralino annihilation in the 
Galactic halo. As opposed to the gamma-rays, positrons 
lose their energy through the propagation in interstel- 
lar space due to the inverse Compton processes and syn- 
chrotron radiation as bended by Galactic magnetic fields. 
For this reason, high energy positrons can only come from 
the region within a few kpc around the Earth. Thus 
positron flux is insensitive to the density profile of the 
dark matter halo and implications of its detection on dark 
matter models are promising [6l|, [63, [63| . 

Propagations of positrons are described by the follow- 
ing diffusion equation, 



d_ 
di 



f(E) = K(E)V 2 f(E) 



_d_ 

~dE 



[b(E)f(E)]+Q(E), 



(40) 



where f(E) denotes the positron number density per unit 
energy, E denotes the energy of the positron in units of 
GeV, K(E) is the diffusion constant, b(E) is the energy 
loss rate and Q(E) is the source term coming from neu- 
tralino annihilation, given as 



Q(E,r) = n 2 (r)(av) 



dE' 



(41) 



where dip/dE denotes the spectrum of the positron from 
single annihilation. The diffusion constant and energy 
loss rate are given by [6I| 



K(E) = 3 x I0 27 [3 06 
b(E) = 10 



E°- e ] cm 2 s-\ 



W E 2 s-\ 



(42) 
(43) 



with E measured in units of GeV. We are interested in 
steady state solution, that is, the solution when l.h.s. of 
Eq. (|4"D|) is set to zero. After solving the diffusion equa- 
tion, the positron flux is given by & e +(E) = (c/4x)f(E) 
with the speed of light c. 
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In order to investigate the detection possibility, we de- 
fine the positron fraction R, 



Re+(E) = 



$e+(E) 



<$> e -{E) + $ e+ {Ey 



(44) 



Here $ e -( e +) includes background flux coming from the 
cosmic ray processes. We use the fitting formula for 
the background positron and electron flux obtained in 
Ref. M,Mi, as 



(prim) 



$ 



$ ( i ec) (£) = 



(E) 



0.16E 



-1.1 



^(sec) 



&r-\E) 



1 + IIE°-V + '3.2E 2 - 15 ' 
0.70E - 7 
1 + 110E 1 - 5 + 600E 2 - 9 + 580£ 42 '' 
4.5£ 07 



1 + 650.E 2 - 3 + 1500S 4 - 



(45) 
(46) 
(47) 



in units of GeV cm s sr . Superscripts (prim) and 
(sec) correspond to the primary and secondary origins of 
them. The reason to use the positron fraction instead 
of the positron flux itself is that the effect of solar mod- 
ulation, which is important for the low energy positron 
flux, is removed by taking the ratio. Another important 
factor comes from the possible local inhomogeneity for 
the dark matter distribution in the Galactic halo, char- 
acterized by the boost factor (BF) , which determines the 
overall normalization of the positron flux [65] . Here we 
conservatively assume BF = 1, which means that the 
homogeneous distribution of the dark matter is assumed. 

Fig. [8] shows typical positron flux for some model pa- 
rameters of MAMSB and mMSB. It is seen that for the 
case of MAMSB, the peak signature will be observed. 
This is because the wino-like neutralino mainly annihi- 
lates into TF-boson pair and they subsequently decay into 
e + , which carry roughly half energy of the primary W + 
boson. In the case of mMSB where the LSP is mostly 
bino-like, positrons are produced through the cascade de- 
cay of the primary annihilation products (ti) , and hence 
the positron flux accumulates at rather low energy re- 
gion. An interesting point is that these dark matters 
may explain the anomaly reported by the HEAT experi- 
ment |66j without introducing any clumpy distribution of 
the dark matter. To explain this anomaly a large annihi- 
lation cross section of dark matter is required in general, 
and the non-thermally produced dark matter naturally 
satisfy it. 

The upcoming experiments such as PAMELA and 
AMS-02 have good sensitivities for a positron energy 
range 10 GeV < E < 270 GeV. We have performed a 
X 2 analysis for investigating detection possibility in these 
experiments following the method of Ref. [63| . The % 2 is 
defined as 



X 



E 



(N° hs - Af G ) 2 



N' 



obs 



(48) 



where JV° bs and Nf G are the number of positron events 
expected background events in the i-th energy bin, re- 
spectively. We chose 22 energy bins as Alog_E = 0.06 
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FIG. 8: Typical positron flux (R(E)) as a function of positron 
energy .E(GeV). Solid line corresponds to MAMSB with 
m 3/2 — 80 TeV and tan/3=10, dot-dashed line corresponds 
to mMSB with M = 450 GeV and tan/3=13. Dotted line 
represents background events. 



for E < 40 GeV and AlogE = 0.066 for E > 40 GeV 
and assumed one year operation. Figs. [9l and 1 101 show the 
resulting x 2 for PAMELA. The \ 2 for AMS-02 becomes 
25 times larger than those for PAMELA. Since 95 and 
99 percent confidence level correspond to x 2 = 34 and 
40 respectively, PAMELA will surely detect dark matter- 
originated positron fluxes in broad parameter region with 
T x < 1 GeV. 
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FIG. 9: Detection posibillity of the positron signal by 
PAMELA in the MAMSB model. We take BF=1, and 95 
and 99 percent confidence level correspond to X — 34, shown 
by dashed line, and 40. The definition of each colors is same 
as Fig. HI 
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FIG. 10: Detection posibillity of the positron signal by 
PAMELA in the mMSB model. We take BF=1, and 95 and 
99 percent confidence level correspond to % 2 — 34, shown by 
dashed line, and 40. The definition of each colors is same as 
Fig. E] 
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FIG. 12: Neutrino flux from the sun in the mMSB model. 
The upper line shows the current upper bound from Super- 
K and the lower one shows the sensitivity of IceCube. The 
definition of each colors is same as Fig. [5] 
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FIG. 11: Neutrino flux from the sun in the AMSB model. The 
upper line shows the current upper bound from Super-K and 
the lower one shows the sensitivity of IceCube. The definition 
of each colors is same as Fig. [5] 



3. Neutrino-induced muon flux from the Sun 

Dark matter particles in the halo scatter off the nucleus 
in the Sun, and then they are trapped and accumulate 
at the center through the lifetime of the Sun. Trapped 
neutralinos annihilate with the enhanced rate because 
of large number density of the neutralinos in the Sun 
[671 168| . The neutralino number accumulated in the Sun 
(N) evolves in a balance between trapped rate and the 
annihilation rate, as 



where Aq = (av)/V denotes the annihilation rate with 
the volume of the Sun V and C© is capture rate by the 
Sun, calculated as [l|, [69| 



C m -3.4 x 10 20 sec _1 




where /9 XS olar and w XS oiar are local density and velocity 
of the dark matter around the solar system, cth and 
one denote the scattering cross section of the neutralino 
with hydrogen and helium respectively, and superscript 
SD(SI) denote spin-(in)dependent component of them. 
These scattering cross sections are limited by the direct 
detection experiments. 

Eq. (|49| is easily solved analytically and we obtain the 
annihilation rate T as a function of time, 



r = \A Q N 2 = ic tanh 2 



VA G C Q t 



(51) 



N = C C ,-A P) N 2 , 



(49) 



Thus we can see that for -y/ ' AqCqI 3> 1, which is valid 
for the case of the Sun (isun — 4.5 Gyr), the annihilation 
rate is simply given by Y — Cq/2 and hence indepen- 
dent of the annihilation cross section of the neutralino. 
Rather, the annihilation rate is determined by the scat- 
tering cross section with nucleus. This is because the 
number density accumulated in the Sun is saturated by 
the balance between the capture rate and the annihila- 
tion rate, and hence the latter is related to the former. 

Once neutralinos annihilate in the Sun, energetic neu- 
trinos produced by the subsequent decay of annihila- 
tion products escape the Sun and reach to the Earth. 
They may be observed at the neutrino detectors such as 
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AMANDA and IceCube, which search muon signals re- 
sulted from high energy neutrino-nucleus interactions in 
the Earth. 

In Fig. II II and 1121 the expected neutrino- induced muon 
fluxes are shown. Here we have taken the threshold 
energy Eth = 1 GcV. The muon flux should be re- 
stricted below 10 3 km~ yr _1 by the current experimen- 
tal bounds, but almost all parameter region is free from 
this constraint. The expected sensitivity of IceCube is 
around $ M ~ 10 2 km _2 yr _1 and only the Higgsino-like 
neutralino in the mMSB model would be promising to 
be detected. The difference from the 7-ray and positron 
signals described in the previous subsections is that the 
high-energy neutrino flux from the Sun is not determined 
by the annihilation cross section of the neutralino, but by 
the scattering cross section of the neutralino with nucle- 
ons. Since the scattering cross section is in general not 
correlated with relic abundance, neutrino-induced muon 
signals are not so enhanced even if non-thermal dark mat- 
ter scenario is assumed. 



IV. CONCLUSIONS 

We have discussed direct and indirect detection signa- 
tures of neutralino dark matters produced non-thermally 
with a very low reheating temperature Td ~ C(l)MeV — 
C(l)GeV. In this scenario, the self-annihilation cross 
section of dark matter should be large enough to account 
for the present relic abundance, within the rage consis- 
tent with bounds provided in [7C|. In SUSY models, such 
a large annihilation cross section is naturally realized for 
the neutralino dark matter with significant wino or Hig- 
gsino components. In the case of bino-like neutralino, 
such a large annihilation cross section can be obtained 
by s-channel Higgs resonance. In both cases, the large 
annihilation cross section leads to the enhancement of 
gamma-ray signals and the positron flux from the dark 
matter annihilation, and it becomes promising to detect 
the non-thermally produced dark matter with T x ^5 lGeV 
by these indirect detection experiments. In other words, 



the indirect detection experiments may give us clues to 
explore the history of the universe with the temperature 
up to 1 GeV. 

Obviously, the consideration with other observations 
is important and essential to make definitive conclusion 
about the non-thermal production scenario. For exam- 
ple, it is known that the large annihilation cross sec- 
tion of the LSP affects the Big-Bang-Nucleosynthesis 
(BBN), and non-thermally produced dark matter with 
T x ^S O(100)MeV may be severely constrained by the ob- 
servation of 6 Li abundance [7l|. And also the combina- 
tion with collider experiments may be most important. 
The upcoming Large Hadron Collider (LHC) experiments 
are expected to discover new particles relevant to the 
EWSB in the standard model, and there may appear 
dark matter candidates. Once such a dark matter can- 
didate is discovered, we may have insight on its produc- 
tion mechanism in the universe by comparing the theo- 
retical calculation of the cross section with cosmological 
and astrophysical observations, such as the dark matter 
abundance and its direct and indirect detection signa- 
tures. As explained in this paper, large indirect detection 
signals are characteristic features for the non-thermally 
produced dark matter, and combined with LHC experi- 
ments we may probe the nature of dark matter by these 
experiments. 

Note added: While finalizing this manuscript, Ref. [72| 
was submitted to the preprint server, which studied simi- 
lar subject to our present work. While they focus on non- 
thermally produced wino and higgsino like dark matter, 
we have also studied bino-like one in the s-channel res- 
onance region, and performed detailed parameter analy- 
ses in MAMSB and mMSB models. The main conclusion 
seems to be consistent with ours. 
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